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Abstract. Halo initial mass functions (IMFs) heavily biased toward 
white dwarf (WD) precursors (~ 1 — 8M Q ) have been suggested as a 
suitable mechanism for explaining MACHO events. However, by simple 
chemical evolution argument, Gibson & Mould (1997; GM97) pointed out 
that such WD-heavy IMFs (wdlMF) cause the overproduction of carbon 
and nitrogen. We re-examine this problem using numerical simulations. 



1. The Code & the Model 

Numerical simulations are performed using an updated version of Kawata (2001). 
Metallicity dependant cooling is incorporated. Star formation occurs when p > 
Pcrit = 2x 10~ 25 gcm~ 3 Sz A-Vi < and follows the Schmidt relation SFR cx p . 
We use two IMFs, the canonical IMF of Salpeter (1955) (sIMF) and wdlMF 
of Chabrier (1996). Stars > 1OM explode as Type II supernova. We use 
metallicity dependent yields from Woosley & Weaver (1995) for Type II SN 
ejecta, which are instantaeously fed back into the ISM. For stars 1 — 8M we 
use the yields of van den Hoek & Groenwegen (1997) and recycling has a time 
delay of 1 Gyr to account for stellar lifetimes. 

We calculate a semi-cosmological disk galaxy formation scenario after Katz 
&; Gunn (1991). Our seed galaxy is an isolated sphere upon which we superim- 
pose small-scale density fluctuations corresponding to a CDM power spectrum. 
Longer wavelength fluctuations are incorporated by enhanced sphere density 
and initial solid-body rotation. We use 9171 gas and dark matter particles for 
a total galactic mass of 5 x 10 11 Mq and examine the chemical evolution of our 
simulated galaxies using the sIMF, wdlMF and a variable IMF (vIMF). For the 
vIMF we use the wdlMF in low metallicity (Z < 0.05^0 ) star forming regions 
whilst we use the sIMF in high metallicity star forming regions. 



2. Results & Conclusions 

Our simulations show that sub clumps of low metallicity stars become read- 
ily apparent by z=2, followed by a continued hierarchical collapse, while high 
metallicity stars are formed primarily at the time of disk formation. These large 
merger events occur in the epoch between z=l.l and z=0.8. High metallicity 
stars preferentially end up in the disk and bulge regions, consistent with White 
& Springel (2000). Table 1 shows that the wdlMF scenario produced a factor 
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of 5 more WDs than the sIMF model. The number of WDs using the vIMF is 
similar to that of the wdlMF in the halo region and midway between the sIMF 
and wdlMF in the disk region. The vIMF yields a local halo WD mass fraction 
of ~ 2% consistent with the observations of Oppenheimer et al. (2001). 

We examine [C/O] vs [Z] for our simulated stars. Observational constraints 
of [C/O] for halo stars ([Z] < —1.5) come from Timmes et al. (1995). Disk 
stars {[Z] ~ 0) have [C/O] ~ 0. Abundances of [C/O] for sIMF correspond 
to observations of both halo and disk stars. For wdlMF [C/O] increases and 
becomes too high once intermediate stars die and their high C yields are recycled. 
vIMF traces wdlMF in the low metallicity region. As metalllicity increases, 
carbon is not overproduced in vIMF nearly as dramatically as in wdlMF. [C/O] 
approaches that for sIMF as [Z] increases further, finally recovering [C/0]~ 
for disk stars. 

The difficulty with reconciling a heavily WD populated halo with abun- 
dances of carbon is consistent with GM97. We find that the overabundance of 
carbon was mitigated by use of a simple variable IMF model in which sIMF 
was used in low metallicity star forming regions and wdlMF was used in high 
metallicity star forming regions. In addition our dynamical simulation reveals 
that vIMF leads to a high population of WDs in the halo in the process of 
hierarchical clustering. 

Table 1. Number densities [pc -3 ] of white dwarves in the 3 models. 
vIMF sIMF wdlMF 

Halo density 7.2 x 10~ 5 2.2 x 10~ 5 9.8 x 10~ 5 
Disk density 2.7 x 10~ 3 1.0 x 10~ 3 6.0 x 10~ 3 
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